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• Abstract 

The instabihty of the ^^O nucleus with respect to decay through the two- 
neutron channel is investigated. It is shown that this isotope, unobserved in 
I the fragmentation experiments, can exist as a narrow resonance in the system 

■ ^^O + 2n. A role of deformation in formation of the neutron-drip line in the 

region A^ ~ 20 is discussed. 
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■ The structure of nucelei near the neutron-drip line possesses interesting features. 

These are: 1) two- neutron instability of nuclei which are stable with respect to 
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O ! one-neutron emission, 2) arising a new region of deformation when approaching 
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the neutron-drip line, 3) change of the concept of nuclear shells at the boundary of 
neutron stability, 4) formation of the neutron halo, 5) clasterization in very neutron- 
rich nuclei. In this article, we pay attention to nucleon stability of oxygen nuclei near 
the neutron-drip line. This region, close to double magic area, has been attracting 
^ ! an attention of nuclear experimentalists an theorists for last decade. 

^ ' The instability of the ^^O nucleus has been observed for the first time at the LISE 

spectrometer at GANIL in the fragmentation of the '^^Ca beam with an energy of 44 
^ ■ Mev per nucleon on a Ta target fl]]. The work p confirmed this result. No events 

. associated with the isotope ^^O were also observed in the recent analysis of oxygen 

fragments by the separator RIPS in the projectile fragmentation experiments with 
using 94.1 MeV per nucleon ^''Ar beam at the RIKEN Theoretical study of 
properties of ^^O was done in within the framework of the self-consistent theory 
of finite Fermi systems [^, H. Estimation of one-neutron separation energy yielded 
Sn — 0.7 MeV, i.e. this nucleus was found to be stable to one-neutron emission. 
This clears the way for the pure two-neutron instability of the ^^O nucleus. Before 
going over to investigation of this subject we concern briefly neighbouring even-even 
oxygen isotopes. Calculations made in for the nucleus ^^O observed in all the 
experiments yielded Sn = 3.6 MeV and = 5.9 MeV. At the same time those 
calculations showed In- and 2n-instability of the nucleus unobserved in ^ ^. 

The two-neutron instability corresponding to positive one-neutron separation 
energy Sn and negative two-neutron separation energy 5*2^ is possible owing to the 
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positive energy of neutron pairing [^]. This energy is usually characterized by the 
even-odd difference A„ of the binding energies of three neighboring isotopes, 
which is given by 
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EhiN+2,Z)-2EbiN+l,Z) + Ef,iN,Z) (1) 



and is positive for even-even nucleus (A^, Z). The obvious relation S2n — '^Sn = —2 An 
follows from eq. (0). This relation indicates that two- neutron instability is possible 
provided that the condition 0<Sn<An is satisfied. Bearing in mind that in this 
region of nuclei A^ ~ 2 ^ 3 MeV one can see that this condition holds for the ^^O 
nucleus. 

To investigate the stability of the ^^O nucleus, we consider the ^"^O -|- 2n system. 
The energies Eg and widths Tg of quasistationary states of this system are deter- 
mined by the poles Eg — iV /2 of the Green's function of two neutrons in the mean 
field of the ^^O nucleus 0. The equation for the two-neutron Green's function can 
be reduced to the equation for the two-neutron wave function \E'(1, 2): 

i7o(l) + H,{2) + 2)] vl/(l, 2) = E^il, 2), (2) 

where the arguments 1 and 2 stand for the sets of neutron spatial and spin co- 
ordinates, Hq is the single-particle Hamiltonian, and V is the effective interaction 
in the two-particle channel. Recent microscopic calculation of this interaction for 
semi-infinite nuclear matter showed its surface nature. The phenomenological 
analysis of even-odd effects in tin isotopic chain showed that density dependence 
of the effective pairing interaction is quite sophisticated. However the parameters 
controlling the shape and the strength of this interaction are fitted for rather heavy 
nuclei not very far from the valley of /5-stability and can hardly be directly used 
in the problem under discussion. That is why we did not claim to calculate the 
wave function ^^(1,2) and the energy E. The aim is to estimate the width of the 
two-neutron state by the order of value. 

To solve eq. (|^), we expand the function 2) in the basis formed by the prod- 
ucts of two eigenfunctions of the Hamiltonian Hq with the total angular momentum 
/=0 and its projection M=0. This expansion has the form 
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vl/(l,2) = j j de,de2C''{e^,e2){reSl)-re,m]'\ (3) 



where v is the set of angular quantum numbers. 

For low energies, the e-dependence of the functions iJsiX) can be separated from 



the coordinate and spin dependences as follows ||T0|, |TT|: 



^,(1) = v/a(^^o(1). (4) 
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According to the theoretical scheme of neutron single-partcle levels of the O 
nucleus, the 2si/2 state is the shallowest bound state, and the 1^3/2 state is the 
lowest quasistationary state. At low energies, the factor 
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is resonantly enhanced for the continuum states because of proximity to the 1^3/2 
quasistationary state with energy e'^ = eQ—i'y/2 in the ^^O nucleus Owing to 
existence of the weakly bound 2si/2 state, the wave functions of the Si/2 continuum 
states are enhanced at low energies by the factor ||TT | 
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where (3=2 /niR'^ and R is of the order of the ^^O radius. For this reasons, we restrict 
the basis in the expansion @ of the wave function \E'(1, 2) to the si/2 and c/3/2 states 
of single-particle continuum. 

To derive the required relations, we follow the method used in [|l^ to obtain the 
relations of the theory of two-proton radioactivity. Substituting expansion with 
the functions V'e(l) factorized in form @) into eq. (H), we find that the energy of the 
quasistationary state satisfies the algebraic equation 
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where 
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9xu = -j dld2V{l,2) {^o'(l) ■^o(2)r{^o(l) -^0(2)}°°, (9) 

and A, z/ = s,d. Substituting expresiions and (^) for the quantity A^{e) into 
eq. (I) and calculating the integrals by rotating the contour to the negative imagi- 
nary axis, we obtain 



M'{E) = 2m 



+ i 



(10) 



M'^{E) 



27^ 



E 



2ei 



-E 



Ti{2ei-EY 



-E 



+ 



25^- 



-E 



In 



-E 



Eq. (0) with M'{E) and M'^{E) specified by eqs. (|ToD and (|TT|) is an algebraic 
equation with complex coefficients. Numerical analysis of eq. (|^ shows that it has 
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a complex root E = Eq — iT /2 in a. wide region of input parameters. While the 
position of this root depends on the values of the parameters, it is located around 
£^0 ~ 0.1 MeV, r ~ 10~^MeV. We reserve a detailed study of this dependence 
for the future publication. The width F ~ 10"'^ MeV of this quasistationary state 
corresponds to a lifetime r ~ h/T ~ 10^^^ s. Since r ^ Tnud ~ 10~^^s the emitted 
neutron pair should be a weakly bound di-neutron state. It can be observed in 
correlation experiments. The similar situation was discussed earlier in connection 



with an analyzis of the /5-delayed multi- neutron emission |T2|, |T3[. However, in that 
case the cascade (n+n) neutron emission strongly dominates and the di-neutron 
(^n) channel is suppressed (e.g. for ^^Na, P2n/ Pn+n < 0.19 [|13|) due to a decay of 
the excited state. In the case of ^^O, the neutron pair is emitted from the ground 
state, so that the situation is "more clear" . 

The spherical basis was used in calculations for neutron-rich oxygen isotopes. 
However, account of deformation (even small, with /?2 < 0.2) could strongly change 
the picture. Indeed, splitting of the single-particle 1^3/2 neutron state results in de- 
stroying the shell = 20 for nuclei near the neutron-drip line [|I^]. Such a sensitivity 
to small deformations should essentially complicate description of neutron-rich nu- 
clei. However, the oxygen isotopes evidently should not be deformed. Otherwise, 
being deformed with P2 > 0.1, the isotope ^^O could obtain additional stability due 
to increase of the binding energy by ~ 1 MeV and should be observed experimentally, 
but this is not the case. 

This is the deformation that seems to give an explanation of existence of the 
isotope ^^F observed in the experiment The calculations |^ of one- and two- 
neutron separation energies for this nucleus based on the self-consistent finite Fermi 
system theory in the sphericl geometry yielded nucleon instability of ^^F. At the 
same time the scenario of sharp shift of the neutron-drip line owing to onset of a 
deformation was suggested several years ago for explanation of the nucleon stability 
of heavy sodium isotopes [|I^. Nuclei with less than half occupation of the level /7/2 
being slightly unbound in spherical calculation {Sn ~ 0) can become stable due to 
lowering the energy of states with asymptotic quantum numbers | [330] and | [321] 
at deformation |T^. The analogous scenario seems to take place for fluorine isotopes: 
the neutron level fY/2 just starts to be occupied in the nucleus ^^F with = 22. 
Following this scenario one could expect nucleon stability of the isotope ^^F as well. 
The problem of onset of a deformation for nuclei near the neutron-drip line will be 
discussed in a separate article, in particular, in connection with an opportunity of 
clasterization in weakly bound neutron-rich systems. 

Authors are grateful to B. V. Danilin, D. Guillemaud-Mueller and M. V. Zhukov 
for valuable discusiions. 
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